Condensation in disordered lasers: theory, 3D+1 simulations and experiments 
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The complex processes underlying the generation of a coherent-like emission from the multiple- 
scattering of photons and wave-localization in the presence of structural disorder are still mostly 
un-explored. Here we show that a single nonlinear Schroedinger equation, playing the role of the 
Schawlow-Townes law for standard lasers, quantitatively reproduces experimental results and three- 
dimensional time-domain parallel simulations of a colloidal laser system. 

PACS numbers: 



Random lasers (RL) are a rapidly growing field of re- 
search, with implications in soft-matter physics, light lo- 
calization and photonic devices Jj, t2] ■ Since the pioneer- 
ing investigations [1, [1| , different groups reported on ex- 
perimental observations, from paint pigments to human 
tissue j5i, 6, 7, 8, 9]. In all of these cases a coherent-like 
narrow spectral line emerges from the fluorescence as the 
pump energy is increased and, in some instances, several 
spectral peaks have been reported d, [l^ . 
In standard single-mode lasers, without structural dis- 
order, the emission linewidth is linked to the electro- 
magnetic energy stored in the cavity by the so-called 
Schawlow-Townes (ST) law [ll|, [T2| • An equivalent law 
for RL is missing. Nevertheless various issues (like the 
statistical properties and the link with spin-glass theory 
0[il[Il[lllll[i3)> were theoretically analysed, while 
the leading model (quantitatively compared with experi- 
ments) is that based on the light-diffusion approximation 
[H, [l^, [l^j mi , which however overlooks the ondulatory 
character of the involved photons. Within a different per- 
spective, RL is due to several localized electromagnetic 
(EM) states put into oscillations in a disordered environ- 
ment (as, e.g., in @, [Hi ) • this framework, it 
is expected that the number of involved modes increases 
with the pump energy and, correspondingly, the spec- 
trum widens. However, exactly the opposite happens and 
this is also accompanied by the shortening of the emitted 
pulse [13, m, HI] . In addition, the fact that strong (or 
Anderson) localization of light sustains the RL action is 
still debated. Ab-initio computational studies were lim- 
ited to ID and 2D geometries [13, 113 j not accounting for 
the critical character of three-dimensional (3D) localiza- 
tion [291. Monte-Carlo simulations neglect interference 
effects[3S[3l[33- 

Here we report on an original theoretical formulation; we 
quantitatively compare its predictions with experiments 
and with the first ever reported 3D-f 1 ab-initio Maxwell- 
Bloch simulations. We show that the RL linewidth is 
ruled by a nonlinear differential-equation, which is the 
equivalent of the ST-law, and is formally identical to the 
nonlinear Schroedinger, or Gross-Pitaevskii (GP), equa- 
tion governing ultra-cold atoms There is hence a 



strict connection between photons in RL and ultra-cold 
bosons; the spectral narrowing observed in RL is thus 
ascribed to a condensation process (3^ of the involved 
electromagnetic resonances. 

Simulations — We consider a vectorial formulation of 
the Maxwell-Bloch (MB) equations [H, [sl]. 21 nonlin- 
ear stochastic partial differential equations are solved by 
finite-difference time-domain (FDTD) discretization on a 
grid distributed on (typically) 256 processors. We model 
an active medium that is infiltrated in the voids of a 
granular distribution of particles obtained by molecular 
dynamics [s^l- We consider 8000 Ti02 particles (aver- 
age diameter 300 nm) with refractive index 2.9 (Fig. [T]). 
The gain bandwidth is 230 nm (^ 1/tg with tg the life- 
time) and the central wavelength is Aq = 590 nm [ujq is 
the angular frequency). Amplification is only present in 
the interstices between colloidal spheres with pump-rate 
(varied by the atomic inversion density Na, see l36j ) 
constant over the ~ 3 ps simulation. The lasing action 
is self-starting from the noise due spontaneous emission 
(SE) modeled as a stochastic term. 

In absence of light-amplification, the response to a single 
cycle pulse 1 fs) at wavelength A — 532 nm, (Fig. 
[1} unveils several spectral peaks corresponding to long- 
living modes. Field spatial distribution (inset of FiglT]) 
is determined by continuous-wave (CW) excitation. 
Then we simulate the RL action: when increasing the 
pumping a coherent field is built from noise. Fig [2^ shows 
the snapshot of the EM energy density in the sample 
middle section. In Fig[2b,c we display the spectra for 
two pumping levels in quantitative agreement with the 
experiments below. The 3D-RL action is mediated by 
several modes with overlapping resonances. In real- world 
samples the RL volume and the number of modes N is 
much larger than that found in our simulations; the out- 
come is a smoother emission profile (Fig|4]). 
Gross-Pitaevskii equation — The RL frequency content 
is given by A{uj), such that |A(w)p is the energy stored in 
the disordered cavity at relative angular frequency lu with 
respect to ujq. The loss coefficient is a{uj), and the gain 
g depends on the whole shape ^(i^), as due to the non- 
linear susceptibility of the resonant medium [l3| • In the 
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FIG. 1: (Color online) Electromagnetic spectrum as obtained 
after a wide-band excitation; the left inset shows the energy 
density in the middle section (A — 586.5 nm, asterisk); a 
sketch of the system is also shown on the right. 
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FIG. 2: (Color online) 3D+1 Maxwell-Bloch simulation of RL; 
(a) energy distribution in the sample middle section (A'^a x 
lO^'^m-^); (b) spectrum for iV^ = 2 x lO^^m"^; (c) as in (b) 
for iVa = 5 X lO^^m-^ 

frequency domain, the oscillation condition "gain=loss" 
reads as (/[^(a;)] = a{u!)A{uj). Limiting for the moment 
to the loss profile a{uj), and following the previous nu- 
merical analysis, one has that, in a disordered system 
sustaining various resonances, Q!(a;) is a smooth function 
interleaved by narrow spectral dips, corresponding to N 
localized (high Q-factor) long- living modes (see figure[T]). 
In this case the loss function a{uj) can be modeled as 



a{uj) — ao(^) 



N 

E 



(1) 



where ao{uj) is the non-resonant smooth loss profile and 
ajiu — LOj) is a sharply peaked (centered at ujj) lineshape 
corresponding to the locaHzed mode j [aj{Lo) is centered 
at w = for later convenience]. Since the RL spectral 
line is limited, we take ao{uj) = oq. For large systems, we 
also expect a huge number of modes 3 with comparable 
properties, we will hence take aj{uj) = aQi,g(w), where 



ctavg is an average resonant line-shape. Therefore the 
oscillation condition becomes 



.g[A(ti;)] = ao{uj)A{Lu) — aj{uj — LUj)A{uj) 

3 

= aoA{uj) - ^ ttavgiuJ - UJj)A{LUj), 



(2) 



where we exploited the fact than aj{u! — cuj) is much 
narrower than A^uj), and hence it "samples" the emission 
spectrum at tOj. In the continuum limit the right-hand 
side of Eq. ^ becomes 

a{uj)A{uj) ^ aoA{Lj) ~ J aavg{i^ ~ ^)A{^)(M. (3) 

We then consider the left-hand side (amplifying part) 
5[A(tj)] of Eq. (21) and we exploit the passive mode- 
locking laser theory [H, [s^ . For a finite gain bandwidth 
with life-time tg, Eq. ([5]) in the time-domain is 
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[ao-0L(i)]a, (4) 
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where we introduced the Fourier transform a(i) 
A{lo) — T\a\ — {1/2it) J a{t) exp{iujt)dt and 0l — 
J'[c(avg]- III Eq. (|4]) go is the small signal gain and 7s is 
the gain saturation coefficient [sl, [3§| . aavg (w) is narrow 
with respect to the gain bandwidth, hence 0L(i) can be 
expanded around t = 0: 4'L{t) = (ao — ^l) [l — (V^i)^] j 
where aL is the average loss for the high-Q modes 
(ao > ul) and iL is their average lifetime. Equation 
(jH) is then cast in a dimcnsionless form a — aofir) 
and t — Tto with Oq — tg^ao — ctL/ls^/gatL and = 

_ ^ 



+ r'^ip + \(p\'^(p = E(p. 



(5) 



The "nonlinear eigenvalue" E is given by 



E ■ 
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(6) 



where -p — gojoih is proportional to the pump energy, 
and K = {tg/tL){ao/aL - 1)^^^. 

It is known that wave-resonances in random systems dis- 
play a distribution of decay times that is bell-shaped 
around some value and with comparable width [40|. 
Equation ([5]) is the oscillation condition for these modes 
with different r (which corresponds to the shift from i^), 
including gain saturation, finite gain bandwidth and the 
mode-coupling due to the overlapping resonances. Equa- 
tion ([5]) [or Eq.(|3])] is identical to the bound state GP 
equation for the ID Bose-Einstcin condensation with an 
external potential ipLit). This shows that a spectral re- 
gion of high Q-modes acts as a trapping potential for 
the energy levels of the excited photons. Frequencies 
tend to be concentrated in this spectral range, as Bose- 
condensed atoms tend to be localized by the external 
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trap [33| . Equation ([5|) displays bell-shaped solutions for 
E > 1 (see, e.g., [41]), this implies a pumping threshold 
for the laser action; the corresponding dimensionless gain 
-2/2- 



pth = 1 + reV2 + ka/4+~^/2 is given by = 1 {pth = 1 
as K « 1). As E>1 {p> Pth,) 
tion of Eq. ® is = 2i/4V£^ - 1 exp( 
spectrum at the threshold (i.e. for p ^ pth) is hence 



an approximated solu- 
V2). The RL 



-{E~ l)exp 



Siuj) = \Aicu)\' = ^ 
with the spectral waist (in frequency lu /2Tr 



2TTtgWth = 




(7) 



(8) 



Eq. ([8]) implies that the RL linewidth at threshold 
is a fraction of the gain bandwidth (~ given by 

^/K,/2^/27r « 1. For a gain bandwidth of 250 nm and 
a RL spike linewidth ^ 0.5 nm (i.e. tijtg = 500) and 
taking for the losses ul ^ ao/lOOO [42], it is k = 0.1 
and 2TrtgWth — 0.2. For E > 1, the spectral profile 
is obtained by the numerical solution of (O; from its 
Fourier transform (fi{ojto) (inset in Fig. [SJs) the normal- 
ized waist Wip[E{p)] and peak p^[E{p)] are determined; 
W = w^^jp^lHg and Sp = t^gPiphs are the corre- 
sponding for S{uj) in real- world units (Figs. [3] and (5]). 
Summarizing, Eq.([5]) relates the pumping p to the non- 
linear eigenvalue i?, which fixes the spectral line-shape 
through Eq.Q Jor Eq.Q]; this equation can be hence 
considered as the equivalent for RL of the ST-law. 
Experimental results — We use a colloidal dispersion of 
Ti02 (Sachtleben Hombitan R611) particles in methanol 
doped by Rhodamine B ( Sigma- Aldrich R6626); the 
packing fraction is 0.2 with average index n = 1.5 (the 
measured mean free path by enhanced back-scattering for 
pure methanol is ^! = 1700 nm at A = 532 nm); the RL 
pump is a 120 ps linearly polarized 10 Hz Nd:Yag laser at 
532 nm and 0.8 mm spot-size. Emission is retrieved by 
a fiber coupled spectrograph (Jobin Yvon, focal length 
140 mm) and a thermoelectrically cooled CCD camera. 
Figure m shows the width (standard deviation) and the 
peak of the spectrum averaged over 100 laser shots versus 
pump energy £; the RL line first narrows and then sta- 
bilizes to a smooth profile. The best-fit with the theory 
(Fig. [5]) furnishes k = 0.14; the threshold pump energy 
{E = 1) is £th — 0.09 mJ. In Fig. [SJa we also display the 
linewidth calculated by a Gaussian fit of the FDTD data 
(the energy axis has been scaled to fit the experiments). 
Conclusions — A theoretical approach based on a non- 
linear bound-state equation, identical to the GP equa- 
tion for BEC, [33] has been shown to quantitatively agree 
with experimentally retrieved laser emission in a colloidal 
dye-doped dispersion of Ti02 particles and with 3D-I-1 
first-principle numerical simulations. RL emission can 
be related to a condensation process of several wave- 
resonances in the presence of disorder; the distribution 
of their decay times playing the role of a temporal trap- 
ping potential. The simultaneous spectral and temporal 
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FIG. 3; (Color online) Theory: Peak spectrum (a) and spec- 
tral waist (b) for various k versus the pumping rate in dimen- 
sionless units. Inset: Spectral profiles for E = 1.1, 10,20. 
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FIG. 4: (Color online) Experimental results: (a) spectra at 
energies 20 /iJ (thin line) and 1000 fiJ, (shaded area is a Gaus- 
sian fit), inset: corresponding enlarged central spectral region; 
(b) unitary-area averaged spectra (100 shots) Vs energy. 
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FIG. 5: (Color online) (a) Left scale, experimentally retrieved 
laser peak spectrum versus pump energy, the line is a best fit 
from theory; right scale, nonlinear eigenvalue (dashed lines 
corresponds to E — 1). (b) Linewidth Vs energy; the contin- 
uous line is the best fit from theory; the FDTD Maxwell-Bloch 
simulations (diamonds) are also shown. 
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narrowing with the number photons in RL is hence cor- 
responding to the spectral and spatial narrowing of the 
BEC wave-function at the condensation. 
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